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ABSTRACT 

WD 1704+481 is a visual binary in which both components are white dwarfs. We 
present spectra of the Ha hne of both stars which show that one component 
(WD 1704+481.2 = Sanduleak B = GR577) is a close binary with two white dwarf 
components. Thus, WD 1704+481 is the first known triple degenerate star. From ra- 
dial velocity measurements of the close binary we find an orbital period of 0.1448d, a 
mass ratio, q — Yf"'''" of <? = 0.70+0.03 and a difference in the gravitational redshifts 

' ^ JW faint ^ ° 

of 11. 5+2. 3km s"^ The masses of the close pair of white dwarfs predicted by the mass 
ratio and gravitational redshift difference combined with theoretical cooling curves are 
O.39+O.O5M0 and 0.56+0.07Mo . WD 1704+481 is therefore also Hkely to be the first 
example of a double degenerate in which the less massive white dwarf is composed of 
helium and the other white dwarf is composed of carbon and oxygen. 
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1 INTRODUCTION 

WD 1704+481 was identified as a pair of white dwarfs of 
similar brightness separated by about Garcscec by Sand- 
uleak and Pesch (1982) from objective prism plates. Spec- 
trophotometry of the pair by Greenstein, Dolez & Vauclair 
(1983) showed that, although their visual magnitudes are al- 
most identical, the SE component (WD 1704+481.1 = San- 
duleak A = GR576, V=14.48, (B-V)=-0.09) is bluer and, 
therefore, hotter than the NW component (WD 1704+481.2 
= Sanduleak B = GR577, V=14.45, (B-V)=0.14). This im- 
mediately suggests that the cooler component must have a 
larger radius which also implies a lower mass for a cool white 
dwarf such as Sanduleak B. Greenstein et al. estimated a 
mass of O.32-O.43M0 for Sanduleak B which is well below 
the typical mass of white dwarfs (O.55M0 , Bergeron, Saffer 
& Liebert 1992). 

Low mass white dwarfs such as Sanduleak B are thought 
to be the result of binary star evolution, in which the evo- 
lution of a star during the red giant phase is interrupted by 
interactions with a nearby star. The physics of this interac- 
tion is complex but it is thought to lead to the stripping of 
the outer hydrogen layers from the red giant in a "common- 
envelope" phase, halting the formation of the degenerate 
helium core and leading to the formation of an anomalously 
low mass white dwarf (Iben & Livio 1993). The hypothe- 
sis that binary star evolution forms low mass white dwarfs 
was confirmed by the discovery by Marsh, Dhillon & Duck 



(1995) of at least 5 short period binary white dwarfs in a 
sample of 7 low mass white dwarfs. 

In this paper we present spectra of the Ha line of Sand- 
uleak B which clearly show that it is a close binary with two 
white dwarf components - a double degenerate star (DD). 
We derive the spectroscopic orbits of both components and 
show that it is likely to be the first known example of a bi- 
nary white dwarf with one white dwarf composed of helium 
and one composed of carbon and oxygen. 



2 OBSERVATIONS AND REDUCTIONS. 

The data for this study come from observations obtained 
with the 2.5m Isaac Newton Telescope (INT) in September 
1998 and the 4.2m William Herschel Telescope (WHT) in 
October 1998. The INT spectra were obtained with the in- 
termediate dispersion spectrograph using the 500mm cam- 
era, a 1200 line/mm grating, a 0.9arcsec slit and a TEK 
charge coupled device (CCD) as a detector at a dispersion 
of 0.39A per pixel. The WHT spectra were obtained using 
the blue arm of the ISIS spectrograph, a 1200 line/mm grat- 
ing, a 0.83arcsec slit and an EEV CCD at a dispersion of 
O.22A per pixel. Integration times varied between 600s and 
1200s. Each observation of the stars was bracketed by ob- 
servations of a copper-neon arc. We set the angle of the slit 
so that spectra of both stars were obtained simultaneously. 
The position angle estimated from this slit angle was 295° 
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for the INT and 290° for the WHT, which agrees well with 
the PA of 289° ± 5° given by Greenstein et al. The separa- 
tion of the stars measured from our INT images is 6.0 arcsec 
with an uncertainty of a few tenths of an arcsecond. The 
seeing was good for all the observations (~ larcsec) so the 
spectra of the stars are clearly separated in all our images. 
Of the 49 spectra, 41 were obtained with the INT. 

The bias level in all the images determined from the 
clipped-mean in the overscan region was subtracted from 
all the images before further processing. For the INT spec- 
tra, several images of a tungsten lamp were combined to 
form a normalized master flat-field image for each night's 
data. Similar images for the WHT spectra show mild fring- 
ing. We therefore combined flat-field images taken immedi- 
ately before and after each observation of the stars in order 
to form a normalized flat-field image for each image of the 
star. Extraction of the spectra from the images was per- 
formed automatically using optimal extraction to maximize 
the signal-to-noise of the resulting spectra (Horne 1986) . Un- 
certainties due to photon statistics are rigorously followed 
through the data reduction process so that reliable uncer- 
tainties are known for every point in the final spectra. The 
arcs associated with each stellar spectrum were extracted 
using the same weighting determined for the stellar image 
to avoid possible systematic errors due to tilted spectra. The 
wavelength scale was determined from a fourth-order poly- 
nomial fit to measured arc-line positions. We calibrated the 
wavelength response of the WHT spectra using observations 
of the standard star BD -1-33° 2642 (Bohlin 1996). To cali- 
brate the wavelength response of the INT spectra we used 
a least-squares fit of a smooth function to the ratio of the 
average WHT spectrum of WD 1740-1-481.2 and the average 
INT spectrum of the same star. The core of the Ha line and 
regions affected by telluric absorption were excluded from 
the fit. We then normalized all the spectra using a linear fit 
to the continuum either side of the Ha line. 

We measured the resolution of the spectra by fitting 
a Gaussian profile by least-squares to a neon arc line at 
6532A in a series of spectra. The full- width at half-maximum 
(FWHM) of the model profile was typically 0.94A for the 
INT spectra and 0.45A for the WHT spectra. The INT spec- 
tra for one night are affected by poor spectrograph focus and 
have a resolution of 1.2A. 



3 MEASUREMENT OF THE RADIAL 
VELOCITIES. 

Our spectra of the cooler component of the visual binary 
clearly show two sharp cores to the Ha line (Fig. |l|) which 
vary in position periodically in a sinusoidal manner, indicat- 
ing the presence of two white dwarfs in a close binary orbit. 
An initial estimate of the period was obtained by measuring 
the radial velocity of the deeper core using a model pro- 
file. We first used a simultaneous least-squares fit to all the 
spectra of four Gaussian profiles with a common position but 
independent depths and widths to determine the model pro- 
file. The shape of the profile was then fixed in a least-squares 
fit to the individual spectra in which only the position of the 
model profile was allowed to vary. These initial velocities are 
affected by the presence of the weaker core, but were good 
enough to determine the orbital period. A periodogram of 



these radial velocities showed a clear peak at 6.91 cycles per 
day, which confirmed our initial impression that the period 
is around 0.145d. 

To measure the radial velocities of the two components 
more precisely we used a simultaneous fit to all the spectra 
of two model profiles, one for each star, in which the position 
of each model profile is predicted from a circular orbit of the 
form 7 -I- K sm{2-n{T — To)/P). In this way we are able to 
determine the shape of the two profiles and the parameters of 
the two circular orbits directly. The (variable) resolution of 
each of the spectra and the effects of smearing due to orbital 
motion are included in the fitting process. There are many 
free-parameters in this fitting process so we used a series of 
least-squares fits in which first the profile shapes were fixed 
while the parameters of the orbit were varied and then vtce 
versa, until we had established values for all the parameters 
which were nearly optimal. Only data within 5000 km s^^ 
of Ha and unaffected by telluric absorption are included in 
the fitting process. We used four Gaussian profiles to model 
the broad wings of the Ha line and the core of star with the 
deeper core and two Gaussian profiles for the other star. A 
polynomial is also included in the fitting process to allow 
for smooth, asymmetric features in the profile. For the final 
least-squares fit the parameters of the profile shapes and the 
orbit were all varied independently. The parameters of this 
final fit are given in Table ^. We designate the star with the 
deeper Ha core component B and the star with the weaker 
Ha core is then component C as we will refer to Sanduleak A 
as component A throughout to avoid confusion. 

We also measured the radial velocities of Sanduleak A 
using the same method employed to measure the initial ra- 
dial velocities of Sanduleak B. There is no significant vari- 
ability in these radial velocities, which have an average value 
of 7A = 0.6 ± 0.3km s'\ 



4 DISCUSSION 

The first obvious result of our analysis is the mass ratio 
q ^ M£. = 0.70 ± 0.03. WD 1704+481.2 is one of sev- 
eral double degenerates (DDs) which have been identified 
among low mass white dwarfs (M < O.49M0 ) where the 
fainter companion is sufficiently young, i.e, hot, for its Ha 
core to be detected and a mass ratio derived (Moran, Marsh 
& Maxted, 1999). These observed mass ratios for DDs are 
shown in Fig. ^ Component C has the weaker core in our 
Ha spectra. The depth of the core does not vary strongly 
with temperature over the temperature range expected for 
components B and C (Koester et al. 1998), so we can con- 
fidently state that component C is the fainter component. 
For comparison, we also plot the results of Han (1998). We 
applied a set of selection criteria to the simulated popula- 
tion of white dwarf binaries to approximate these selection 
effects, namely, that the younger white dwarf is less mas- 
sive than O.49M0 and that the older white dwarf is younger 
than 1 Gyr. The mass ratio distribution for white dwarf bina- 
ries predicted by the models of Iben, Tutukov & Yungelson 
(1997) using similar selection criteria but their own "stan- 
dard model" is quite similar, i.e. it shows a single peak near 
q = 0.7. WD 1704-1-481.2 lies near the peak of the theoreti- 
cal distribution, which is a success for the theoretical model. 
What is less clear is why the majority of DDs with measured 
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Figure 1. Examples of our WD 1704+481.2 spectra and the least-squares fit used to derive the orbital parameters of the binary. The 
spectra (which are in no particular order) are offset vertically from each other for clarity and the fit is shown as a smooth line. The 
regions of the spectrum aff'ected by telluric absorption and excluded from the fit are plotted with faint lines. 
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Figure 2. The predicted mass ratio distribution for DDs with 
measured mass ratios from the model of Han (1998). Measured 
mass ratios for double degenerates, including WD 1704+481.2 
(marked with an asterisk) are indicated by error bars. 



Figure 3. The predicted gravitational redshift difl'erence, Avgr, 
as a function of the component masses (solid diagonal line) and 
its uncertainty due to the uncertainty in the mass ratio (dashed 
diagonal lines). The observed value of Avgr (solid horizontal line) 
and its uncertainty (horizontal dashed line) are also shown. 



mass ratios do not have mass ratios which accord with the 
models. 

We can use the observed difference 7b — 7c = — ff.5± 
2.3km to determine the masses of the stars as foUows. 



The gravitational redshift of the light from the stars is 
given by v^^ = 0.635(M/Mq )/(i?/R0 )kms"^ Since one 
star is more massive than the other and more massive 
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Figure 4. The predicted mass distributions for the brighter 
component (horizontal hatching) and fainter component (diag- 
onal hatching) of DDs with measured mass ratios from the model 
of Han (1998). The observed masses for the brighter component 
(asterisk) and fainter component (diamond) of WD 1704+481.2 
are also shown. Note that the uncertainties shown are correlated, 
i.e., an increased mass for the brighter component implies an in- 
creased mass for the fainter component. 



Table 1. Orbital parameters of the final least-squares fit to 
the spectra of WD 1704-1-481.2. The full-width at half-minimum 
(FWHM) and depths (D) of each of the Gaussian profiles used to 
model the Ha line are also given. Other symbols are defined in 
the text. 



Circular orbit 




HJD(To) — 2451000 


77.4070 ± 0.0002 


Pid) 


0.1447864 ± 0.0000025 




— 13.0 ± 0.8 


7^B(kms~^) 


— 135.4 ± 1.3 




— 1.4 ± 2.2 


KnCkxn s~^'l 


94 5 -1- 3 8 


Component B 




FWHMi (A) 


0.86 ± 0.09 


Di 


0.147 ± 0.009 


FWHM2 (A) 


5.8 ± 0.4 


Da 


0.078 ± 0.004 


FWHM3 (A) 


21. 6± 1.1 


D3 


0.116 ± 0.007 


FWHM4 (A) 


55.0± 1.7 


D4 


0.138 ± 0.005 


Component C 




FWHMi (A) 


2.5 ± 0.2 


Di 


0.063 ± 0.003 


FWHM2 (A) 


16.7 ± 1.1 


D2 


0.055 ± 0.004 



white dwarfs have smaller radii (for a given temperature), 
the more massive component will appear to have a larger 
(more positive) systemic velocity, as is observed. To de- 
termine the radius of component B we use the models of 
Althaus & Benvenuto (1997) for helium white dwarfs. We 
assume a temperature of 9500K for component B (Green- 
stein, Dolez & Vauclair 1983) and use parabolic interpola- 
tion to find the radius. Note that a lOOOK difference in tem- 
perature changes the radius derived by less than one per- 
cent so the temperature assumed has a negligible effect on 
our results. We assume that component C is slightly cooler 
(8500K) and use the mixed-composition models of Wood 
(1990) to determine its radius. In Fig. ^ we show the pre- 
dicted gravitational redshift difference, Augr as a function 
of the mass of each component. The observed value of Augr 
is also indicated and we see that this observed value implies 
Mb = 0.39 ± O.O5M0 and Mc = 0.56 ± O.OTM© . White 
dwarfs more massive than about O.SM© will have passed 
through a helium-burning stage and have a carbon-oxygen 
composition, so these masses suggest that component B is a 
white dwarf composed of helium and component C is com- 
posed of carbon and oxygen. In Fig. ^ we see that the masses 
of the two components are roughly what would be expected 
from the theoretical models of Han (1998). 

The mass of component C is typical for single white 
dwarfs (Bergeron, SafTer & Liebert 1992) and component A 
has the same gravitational redshift. The change in radius 
with temperature in this regime is very small, so we expect 
component A has a normal mass. Component A is ^200 AU 
distant from components B and C (Greenstein et al. 1983) 
and the initial separation of components B and C must have 
been smaller than ~ 1 AU for a common-envelope phase to 
have occurred. The ratio of the orbital periods is sufficiently 
large that we need not worry about dynamical instability of 
the orbits (Kiseleva, Eggleton & Anosova 1994) and so we 
can be confident that the evolution of the inner binary will 



have been unaffected by the presence of component A and 
vice versa. 

From the parameters of the orbit we find A/b sin'^ i = 
0.075 ± O.OO3M0 and Mc sin^' i = 0.108 ± O.OO6M0 which, 
combined with the masses derived above, yields an inclina- 
tion of 61°. The separation of components B and C is only 
O.74R0 but the small radii of white dwarfs (~ O.OIR0 ) rules 
out the possibility of observable eclipses in this binary. 



5 CONCLUSION 

We have shown that WD 1704-f 481 is a hierarchical triple 
star in which all three components are white dwarfs, i.e., 
a triple degenerate star. The outermost star, component A 
appears to be a typical white dwarf. Components B and C 
are a close binary with an orbital period of only 0.145d. 
The mass ratio (q = ^^^^psi^ ^ "a ^ o.70 ± 0.03) for 
the inner binary is the first measured mass ratio for a dou- 
ble degenerate which is close to the peak of the mass ratio 
distribution predicted by theoretical models. Similarly, the 
masses of components B and C derived from the difference 
between their gravitational redshifts (A/b = 0.39 ± O.O5M0 
and Mc ~ 0.56 ± O.OTM© ) appear near the peaks of the 
theoretical mass distributions. These suggest that the com- 
ponent C is a typical white dwarf composed of carbon and 
oxygen and that component B is composed of helium. 
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